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p ABSTRACT 

O ' The NIR Call triplet absorption lines have proven to be an important tool for quantitative spectroscopy of individual red giant 

\ branch stars in the Local Group, providing a better understanding of metallicities of stars in the Milky Way and dwarf galaxies and 

C/3 , thereby an opportunity to constrain their chemical evolution processes. An interesting puzzle in this field is the significant lack of 

extremely metal-poor stars, below [Fe/H]=-3, found in classical dwarf galaxies around the Milky Way using this technique. The 
question arises whether these stars are really absent, or if the empirical Call triplet method used to study these systems is biased 
in the low-metallicity regime. Here we present results of synthetic spectral analysis of the Call triplet, that is focused on a better 
understanding of spectroscopic measurements of low-metallicity giant stars. Our results start to deviate strongly from the widely-used 
and linear empirical calibrations at [Fe/H]<-2. We provide a new calibration for Call triplet studies which is valid for-0.5>[Fe/H]>- 
4. We subsequently apply this new calibration to current data sets and suggest that the classical dwarf galaxies are not so devoid of 
extremely low-metallicity stars as was previously thought. 

Key words. Stars:abundances - Galaxies:dwarf - Galaxies:evolution - Galaxies:Local Group - Galaxy:formation 

o . 

1. Introduction and Outline galaxy, as they represent the most pristine (and probably thus 

the oldest) stars in the system. For example, in a scenario where 
To understand galaxy evolution it is critical to understand how t^ere are no low-metalhcity stars, theories invoking some kind of 
. ^ the metalhcities of stars m physically different environments de- pre-enrichment prior to the star formation epoch become more 
velop with time. Low- and high-resolution spectroscopic stud- plausible. In general, detailed chemical abundances of old stars 
H ; les of individual stars in the Milky Way and dwarf galaxies are ^an provide valuable information about chemical evolution pro- 
. ."^ , crucial in this context since they can provide measurements of messes as they are likely to be polluted by only a single or very 
overall metaUicity and detailed abundances of a range of ele- supernova. A comparison of low-metallicity tails in different 

ments (e.g.. Freeman & Bland-Hawthorn 2002, and references galaxies also provides information about the origin and evolution 
therein). Studies of this kind have revealed interesting differ- systems with different masses and different star formation his- 
ences between the stars in the MiUcy Way halo and the dwarf tories 
spheroidal galaxies (dSphs) surrounding it dTolstov et al.ll2009[ 

and references therein). An intriguing population of stars to From relatively time-consuming but insightful high- 
study in different environments are the extremely metal-poor resolution (HR) studies of small samples of bright red giant 
stars ([Fe/H]< -3). The existence, or lack of, these stars reveals branch (RGB) stars in nearby dSphs we know that the chem- 
valuable information about the chemical evolution history of a ical signatures of individual elements in the dSph stars can 



* Using observations collected at the European Organisation for 
Astronomical Research in the Southern Hemisphere, Chile proposal 
171.B-0588. 



be distinct from the stars in the Galaxy (e.g., Shetrone et aL 


2001, 2003; Tolstov et al. 2003; Venn et al. 2004; Fulbriaht et a" 


2004; Kochetal. 


2008dlbl: ICohen & Huand 120091; lAoki et al. 


12009: ,Frebel et al. 


i2010h. So far all high-resolution observa- 
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tions of extremely metal-poor stars in dSphs have been stud- 
ied in high-resolution as part of follow-up studies o n previ- 
ous low-resolution samples (C ohen & Huang 2009; Ao ki et al] 
l2009HFrebel et al.ll2009Ll20ld Tafelmeyer et al. in prep.). These 
stars are still sparse and at the moment mostly found in the 
ultra-faint dwarf galaxies. The high-resolution studies of clas- 
sical dSphs which are not follow-up studies, have first targeted 
the inner galactic regions, which are often more metal-rich, and 
are hence not optimized to specifically target metal-poor stars. 

Additionally, recent low-resolution (LR) studies enable us 
to determine overall metallicity estimates for much larger sam- 
ples of RGB stars in both classical dSphs (e.g., [S untzeff" et aD 
1993'; 'Tolstoy et al."2004; 'Pont et al.'2004; BattagHa et al. 200a 
^unoz et al. 2006; Koch et al. 2007a b; BattagHa et al. 20081 
Shetrone et all f2009t [Walker et al. 2009a; Kirbv et al. 2009), 



ultra-faint galaxies (e.g..iSimon & Geha 2007; Kirbv et alj2008t 
iNorris et alJ lIOOSi: IWalker et al.1 2009b; Koch et al.)|2009D and 
even the more distant and isolated dwarf irregular galaxies (e.g., 
iLeaman et al. 2009). From the larger numbers of stars studied 
in the low-resolution studies for the classical dSphs (typically 
several hundred per galaxy) one would statistically expect to 
find some RGB stars with [Fe/H]< -3, if the distribution of 
metalhcities in these systems follows that of the Galactic halo 
(iHebni et alJ2006) . However, one of the compelling results from 
studies of large samples of RGB stars is a significant lack of 
stars with metallicities [Fe/H]< -3 in the classical dSph galaxies 
Sculptor, Fornax, Carina and Sextans compar ed to the metalUc - 
ity distribution function of the Galactic halo dHelmi et alj|2006h . 
These metallicities are inferred from th e Une strengths of the 
Call NIR triplet (CaT) Unes. Recentlv. iKirbv et al.( (l2009h re- 
ported to have found a RGB star in Sculptor with a [Fe/H] value 
as low as -3.8 using a comparison between spectra and an ex- 
tensive spectral library. Several extremely low-metallicity stars 
have already been discovered in the ultra-fa int dwarf galaxies 
using either this te chnique or other indicators dKirbv et alJ2008t 
iNorris et al.ll2008h . In this paper we investigate whether the lack 
of low-metallicity stars in the classical dwarf galaxies could be 
a bias due to the Ca II NIR triplet (CaT) indicator used to deter- 
mine metallicities from low-resolution spectra. 

The CaT lines have been used in studies over a wide range of 
atmospheric parameters and have been applied to both individ- 
ual stars and integrated stellar populations in different environ- 
ments (see Cenarro et al. 2001, 2002, and references therein). 
In this paper we focus on the use of the CaT lines to deter- 
mine metallicities of individual RGB stars. The CaT region of 
the spectrum has proven to be a powerful tool for metallicity 
estimates of individual stars. The three CaT absorption lines 
(A8498, A8542, and .i8662 A), which can be used to deter- 
mine radial velocities and to trace metallicity (usually taken as 
[Fe/H]), are so broad that they can be measured with sufficient 
accuracy at a moder ate resolution. As was already noted in pi- 



oneering work (e.g., Armandroff & Zinn" 19881 : [Olszewski et alJ 



ll991l;lArmandroff & Da Costa,1991,) . there are numerous addi- 
tional advantages to the use of the CaT lines as metallicity in- 
dicator. For instance, the calcium abundances are expected to 
be largely representative of the primordial abundances of the 
star, since, contrary to many other elements, they are thought to 
be unaltered by nucleosynt hesis processes in intermediate- an d 
low-mass stars on the RGB divans et al.ll2001l:ICole et alj|2004l) . 
Also, the NIR wavelength region is convenient: in that the red gi- 
ants emit more flux in this part of the spectrum than in the blue 
and the spectrum is very flat, which facilitates the definition of 
the continuum level to measure the equivalent widths of the line. 



However, the breadth of the lines also has disadvantages. 
Because the lines are highly saturated their strength, especially 
in the core of the line, depends strongly on the temperature struc- 
ture of the upper layers of the photosphere and chromosphere 
of the star, which means that complicated non-local thermody- 
namic equilibrium (non-LTE) phy sics has to be used to model 
the line correctly (iCole et alJ 20041) . Also, the lines do not pro- 
vide a direct measurement of [Fe/H], although it has been shown 
that [Fe/H] affects the e quivalent width of the lines more than Ca 
dBattagliaet al.l2008bh . The abundance of Ca and other elements 
do still afi^ec t the equivalent width s which will not always trace 
only [Fe/H] dRutledge et al.l[T997bD . 

Early investigations of the CaT concentrated mai nly on 
their sensitivi t y to s urface gravity ( Spinrad & Tavlor 1 19691 
1971; Cohen 1978; Jones et al.1 Il984f) '. It was realized by 
Armandroff & Zinn. (J988) that a more metal-rich RGB star 
should have stronger CaT lines, because it has both a greater 
abundance of Fe (and also Ca) in its atmosphere and a lower 
surface gravity. They empirically proved this relation by mea- 
suring for the CaT lines their integrated equivalent width (EW) 
in several globular clusters with k nown [Fe/H]. Applying this 
method to individual RGB stars, lOlszewski et al.l dl991h no- 
ticed that the metallicity sensitivity of the CaT line index is im- 
proved by plotting it as a function of the absolute magnitude of 
the star. At a fixed absolute magnitude, higher metallicity RGB 
stars will have lower gravity and lower temperatures, which both 
strengthen the CaT lines. A f urther development of the method 
d Armandroff^ & Da CostallT99 1) was to plot the equivalent width 
as a function of the height of the RGB stars above the horizon- 
tal branch (HB) in V-magnitude {V - Vhb)- In this way, the re- 
quirements of a distance scale and a well-determined reddening 
are avoided. They found a linear relation between [Fe/H] and 
a "reduced equivalent width" (W ), which incorporates both the 
equivalent width of the two strongest lines at /18542 and /18662 A 
(also written as EW2 and EW3 in the remaining of this paper) 
and V - Vhb- This enables a direct comparison between RGB 
stars of different luminosities. This method has been extensively 
tested and proven to work on large samples o f individual RGB 
stars in globular clusters (e.g.. lRutledge et al][T997at,b) . 

Additionallv. lCole et all d2004|) showed tiiat the eff'ect of dif- 
ferent ages of RGB stars is a negligible source of error for metal- 
licities derived from the CaT index. This paved the way for 
the use of the CaT metallicity indicator on populations of RGB 
stars which are not coeval, among which are the Local Group 
dwarf galaxies. A direct detailed comparison between the low- 
resolution CaT metallicities and high-resolution measurements 
for large samples of RGB stars in the n earby dwarf galaxies 
Fornax and Sculptor is given by BattagHa et al.l (|2008b). They 
concluded that the CaT - [Fe/H] relation (calibrated on globular 
clusters) can be applied with confidence to RGB stars in com- 
posite stellar populations over the range -2.5 <[Fe/H]< -0.5. 

In this paper we study the behavior of the CaT for [Fe/H] 
< -2.5, comparing stellar atmosphere models with observa- 
tions. We have determined a new calibration including this low- 
metallicity regime. The paper is organized as follows. In Sect. 
|2] we further clarify and physically motivate the need for a new 
calibration at lower metallicities using simple synthetic spectra 
modeling. In Sect.[3]we describe our grid of synthetic spectra we 
use to investigate the CaT lines. We then analyze and calibrate 
this grid in two regimes: -2.0 > [Fe/H] > -0.5 in Sect. |4l and 
[Fe/H] < -2.5 in Sect.|5| In Sect. |6] we present a new calibration 
valid for both metallicity regimes. Additionally, we investigate 
the role of a elements in Sect. |2| and the implications of the new 
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calibration on the results of the Dwarf Abundances and Radial 
velocities Team (DART) survey (iTolstov et al.ll2006h in Sect. [8] 



2. The CaT at low metallicity 

The so-called CaT empirical relation connects a linear combi- 
nation of the equivalent widths of the CaT lines (the exact form 
can vary between different authors) and the absolute luminos- 
ity of the star, often expressed in terms of the height of the star 
above the HB of the system, to its rFe/H1 va lue. The empirical 
relation described in 'Battaglia et al.' ('2008b) in their equations 
16 and 1 1 are given as equations[T]and|2]below. 



[Fe/H] = -2.81(+0.16) + 0.44(+0.04)H^ 



where: 



W = EWi + EWi + 0.64(+0.02)(y - Vhb) 



(1) 



(2) 



At present this linear relation between metallicity and CaT 
equivalent widths is also used to infer metallicities for stars out- 
side the calibrated regime (-2.5 <[Fe/H]< -0.5). However, the 
assumption that the relation continues outside of this regime has 
not been accurately checked. It is very clear that the linear re- 
lation cannot hold down to extremely low metallicities, since at 
a certain point Eqs[T| and |2] infer that the equivalent widths will 
have negative values. Negative equivalent widths for absorption 
lines is obviously physically meaningless. 

It is also expected that the shape of the lines, and therefore 
the relation of equivalent width with metallicity, changes as the 
metallicity drops. The equivalent width of the CaT lines is dom- 
inated at higher metallicities by the wings of the line. At lower 
metallicities these wings weaken and/or even disappear and the 
equivalent width of the line becomes more dominated by its core, 
as illustrated in Fig.[T]from synthetic spectra. Since the core and 
the wings of the lines originate from different parts of the curve 
of growth the behavior of both is unlikely to be described by one 
linear equation. The change in relative strength of the core and 
wings of the hne changes the sensitivity of the line to the Ca 
abundance and to the overall metallicity in general. This under- 
standing of the physical process motivates our re-calibration of 
the relation between CaT equivalent width and [Fe/H] at lower 
metallicities. 



3. A grid of Models 

In order to investigate and quantify the behavior of the CaT lines 
at low metallicity we define a grid of model spectra created us- 
ing the publicly available and recently revised version of the 
(OS)MARCS models (e.g., Gustafsson et al. 2008; Plezi l2008l) 
and the Turbospectrum program (Alvarez & Plez 1998) updated 
consis tently with the latest release of MARCS Gustafsson et alj 
(l2008l) . The model spectra cover a range of effective tempera- 
tures, gravities, metallicities ([Fe/H]) and enhancements of the 
a elements (taken as O, Ne, Mg, Si, S, Ar, Ca, and Ti in the 
MARCS models (iGustafsson et aO2008h ). The models use a ID 
spherical symmetric approach. Our linelist is created using the 
Vien na Atomic Line D atabase (VALD) for the atomic species 
(e.g., iKupka et alJl2000 ). Additionally we also model the con- 
tribution from CN molecular lines (Spite et al. 2005, Plez priv 
comm.) and TiO m oleculer lines in the cooler (Tefj < AQQOK) 
models (lPlezlll998h . 
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Fig. 1. The change of the shape of the CaT line at A 8542 Awith 
changing metallicity (metallicity is decreasing as the line nar- 
rows). Non-LTE effects are not taken into account. These lines 
are obtained from synthetic spectra from MARCS models with 
Teff = 4500 K and log(g) = 1.5 as described in the text. The flat- 
tening in the core at higher metallicities arises because the mod- 
els do not cover small enough optical depths in the outer layers, 
where the core of the line is formed. This is a problem within the 
modeling itself, and not related to the non-LTE effects which are 
described in Sect. 13.21 Because of the broad wings, the area of 
the core of the line that is missing is a negligible fraction of the 
total equivalent width. 



3.1. Parameters 

The parameters [Fe/H], [a/Fe], Teff and log(g) for MARCS 
model atmospheres in our grid are given in Table [1] For the 
upper part of the table all the models are provided on the web 
by the MARCS team. We interpolated to a finer grid, including 
the values given in the lower part of Table [1] using the inter- 
polation tool provided on the MARCS website by T. Masseron. 
All atmospheric models have a mass of \Mq and microturbu- 
lent velocities (^) of 2 km s"'. In calculating the synthetic spec- 
tra with Turbospectrum we use microturbulent velocities which 
vary slightly dependent on the gravity of the atmospheric model 
in accordance with observed variations of microturbulent veloc- 
ities in halo stars. The values we use range from ^ = 2.3 km 
S-' for the tip of the RGB (log(g) = 0.5) to ^ = 1.7 km s ' 
near t he HB (log(g)=2.5) in accordance with the results from 
iBarklem et al.l (l2005h . This slight change in microturbulence can 
significantly alter the equivalent width of the two strongest CaT 
lines in the metal-poor regime and the difference will be largest 
at the tip of the RGB. 



3.2. Non-LTE effects 

Although all model atmospheres and synthetic spectra are cal- 
culated assuming local thermodynamic equilibrium (LTE), we 
take effects of departures from LTE into account, because they 
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Table 1. The parameters for the grid of models used. 



[Fe/H] 


[a/Fe] 


Teff 


log(g) 


MARCS model atmospheres 


-0.25 


+0.0, +0.1 


3800, 3900, 4000, 4250, 4500, 4750, 5000 


0.5, 1.0, 1.5,2.0, 2.5 


-0.50 


+0.0, +0.2 


3800, 3900, 4000, 4250, 4500, 4750, 5000 


0.5, 1.0, 1.5,2.0, 2.5 


-0.75 


+0.0, +0.3 


3800, 3900, 4000, 4250, 4500, 4750, 5000 


0.5, 1.0, 1.5, 2.0, 2.5 


-1.00, -1.50,-2.00 


+0.0, +0.4 


3800, 3900, 4000, 4250, 4500, 4750, 5000 


0.5, 1.0, 1.5, 2.0, 2.5 


-2.50, -3.00 


+0.4 


3800, 3900, 4000, 4250, 4500, 4750, 5000 


0.5, 1.0, 1.5, 2.0, 2.5 


-4.00, -5.00 


+0.4 


3800, 3900, 4000, 4250, 4500, 4750, 5000 


1.0, 1.5,2.0, 2.5 


Interpolated model atmospheres 


-1.25, -1.75 




4125, 4375, 4625, 4875 


0.75, 1.25, 1.75, 2.25 


-2.25, -2.75 




4125, 4375, 4625, 4875 


0.75, 1.25, 1.75, 2.25 


-3.50, -4.50 




4125, 4375, 4625, 4875 


1.25, 1.75, 2.25 



can be significant since the lines are so highly saturated. The ef- 
fec t of departures f r om LT E on these lines was first investigated 
by iJorgensen et al.l (Il992h . but only for [Fe/H] > -1. A more 
extensive study is performed bv lMashonkina et al] (l2007b . They 
note that in the CaT the non-LTE effects are revealed only in the 
Doppler core, which is strengthened. 

Because the broad wings of the CaT lines are decreasing sig- 
nificantly with metallicity, as shown in Fig. [T] the effect of de- 
partures from LTE, which only affect the core, have most im- 
pact on the equivalent width determination at low metallicities. 
Therefore, it is crucial to take non-LTE effects into account in or- 
der to understand the behavior of the CaT lines at low metallic- 
ity. We perform non-LTE calculations using the model atom pre- 
sented by Mashonkina et al. (2007), which contains 63 levels of 
Ca I, 37 levels of Ca II and the ground state of Ca III. Non-LTE 
level populations and synthetic spectra were determine d with re- 
cent versions of the codes DE TAIL and SURFACE (iGiddingsl 
[1981'; 'Butler & GiddingsllMl- We chose ATLAS9 atmospheric 
models (Kurucz 1993) as the input models in the non-LTE com- 
putations. Thus, we computed ATLAS9 model atmospheres ex- 
actly for a given set of stellar parameter s and metallicities us - 
ing a Linux version of the ATLAS9 code (ISbordone et al.ll2004l) . 
and adopting the new Opacity Distribution Functions (ODFs) of 
[Castelli & Kurucz (2003). We also adopted S'h = 0.1 as the scal- 
ing factor of the inelastic collisions with hydrogen atoms in the 
non-LTE computatio ns. For further details on the non-LTE com- 
putations we refer to lMashonkina et al] (l2007h . 

We determined the equivalent widths by integrating normal- 
ized fluxes of the broad Ca II triplet non-LTE profiles at 8542 
and 8662A. Figure |2] shows the ratio of LTE EWs to non- 
LTE EWs from this modeling for both these CaT lines in sep- 
arate panels. At high temperature and [Fe/H]=^.0 the ratio 
goes down in both lines to ~0.7, which means that just ~70% 
of the line is modeled in LTE and non-LTE effects are thus 
very important. We determine a best-fitting relation as a func- 
tion of metallicity and temperature of the model using the IDL 
function MPFIT2DFUN (Markwardt 2009), which performs a 
Levenberg-Marquardt least-squares fit to a 2-D function, in com- 
bination with a statistical F-test to identify the best fit. The best- 
fitting relations obtained separately for the two strongest CaT 
lines are shown as dashed gray lines in Fig. |2] and given in Eqs. 
Hands 



leff 



EW(8542A)lte/EW(8542A)„„„.lte = 0.563 + 0.397 



leff 



4500(K) 



(0.365 - 0-323 ^^J|j^)[Fe/H] - 0.0203[Fe/H]2(3) 



EW(8662A)lte/EW(8662AWlte = 0.692 + 0.292^^^^^^^ 

- (0.301 - 0-288^^J^)[Fe/H] - 0.0163[Fe/H]2(4) 

For our finer grid of synthetic spectra from MARCS model 
atmospheres described in the previous paragraph we use these 
equations to determine the ratio of LTE to non-LTE equivalent 
widths for the two lines in each individual model. By dividing 
the individual equivalent widths in the grid (which are all cal- 
culated in LTE) by the factor for the corresponding line before 
adding the two strongest lines together, we correct each of the 
grid points for non-LTE effects. 

4. The CaT lines at -2.0 < [Fe/H] < -0.5 

4.1. The empirical relation 

The empirical relation between CaT EW, absolute magnitude 
and [Fe/H] is very well studied in globular clusters, i.e., in the 
metallicity range between [Fe/H]a;-2.3 and [Fe/H]ss-0.5. This 
well-known relation can thus be used to test our synthetic spec- 
tra at these metallicities. We make a comparison between the 
CaT lines from atmospheric models between -0.5 < [Fe/H] < 
-2.0 and the best fit linear relat ion for globular cluster s as given 
in Eqs. □ and H] in Sect.l2l(from lBattaglia et"al]|2008bl) . 

Some approximations are needed to enable a CaT analysis 
with synthetic spectra which is comparable to observations. 

1. The synthetic spectra are all degraded to a resolution of 
R=6500, the resolution of VLT/FLAMES used in Medusa mode 
with the GIRAFFE LR (LR8) grating as was used in the ob- 
servational determination of the empirical relations given above. 
The equivalent widths for the two stron gest CaT lines are mea - 
sured using the same fitting routine as in Battaglia et al.l (l2008bl) : 
a Gaussian fit with a correction which comes from a comparison 
with a the summed flux contained in a 15 A wide region centered 
on each line. The correction is necessary to account for the wings 
in the strong lines which are distinctly non-Gaussian in shape. 
Since the CaT lines can have a variety of shapes for a range of 
metallicities, as shown in Fig.[Tl it is in general not advisable to 
fit them using a single profile. The disadvantage of using just nu- 
merical integration of the observed spectra is that there are also 
some weaker lines in this wavelength range that may vary differ- 
ently with changi ng stellar atmospheri c parameters than the CaT 
lines themselves dCarrera et al.ll2007l) . Taking these considera- 
tions into account, we thus use a combination of both methods. 
Weak nearby lines in the spectrum still introduce a small depen- 
dence of the measured CaT EW on resolution though, since they 
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are more likely to be absorbed into the Gaussian fit at a lower 
resolution. This dependence is already present the lowest metal- 
licity in our grid and gets stronger at higher metallicity, due to 
the increasing prominence of the non-Gaussian wings. 

Some of the more prominent weak lines which can be present 
in the wings of the two strongest CaT lines are the hydrogen 
Paschen lines. Their strength mainly depends on the tempera- 
ture of the star (the hotter the stronger) and also to a lesser ex- 
tend on its gravity (increasing with decreasing gravities). In the 
more metal-rich part of our grid the Paschen lines fall within the 
broad wings of the CaT lines and have a direct effect on the EW 
measurement of the CaT line. However, within the range of pa- 
rameters we use in our models, the maximum contribution of the 
Paschen lines to the CaT EW measured is 39 mA, which is negli- 
gible compared to the total CaT EW. In the more metal-poor part 
of the grid the CaT lines are narrower and well separated from 
the Paschen lines. Nonetheless, the Paschen lines can still influ- 
ence the CaT EW by effecting the placement of the continuum. 
Also this effect we find to be negligible at our grid parameters, 
at maximum the CaT EW is changed by 3.5%. 

2. Not all the models in our grid, each of them a particu- 
lar combination of effective temperature, gravity and metallicity, 
will represent real stars on the RGB. To determine which mod- 
els best compare t o real stars we use two s ets of isochrones, the 
BaSTI isochrones ('Pietrinferni et al."2004) and the Yonsei-Yale 
set (e.g.. lYi et al.l l2001: Demarque et al. 2004). Both sets can be 
interpolated to obtain exactly the desired metallicity and age for 
a particular isochrone. The Yonsei-Yale set has the advantage 
that they go down to [Fe/H] = -3.6 at [or/Fe] =-1-0.4, whereas the 
lowest value for the BaSTI set is [Fe/H] - -2.6. It is well known 
that different sets of isochrones do not always give identical re- 
sults. On top of that, different ages give (slightly) different pa- 
rameters for the RGB stars. In Fig.|3]the grid of models in effec- 
tive temperature versus gravity is plotted along with the relations 
from the theoretical isochrones. In the two panels the differences 
due to age, metallicity and the use of different sets of isochrones 
are illustrated. Based on Fig.|3]we decide to linearly interpolate 
in log(g) space between the models that are as close as possible 
to the isochrones and add another 0.25 in log(g) space on each 
side to account for uncertainties within the isochrone models as 
well as the differences between the isochrones of different ages. 
These uncertainties are shown as error bars on the equivalent 
widths from synthetic spectra. Note that these represent max- 



Fig. 2. The ratio of LTE to full non-LTE equiv- 
alent widths for the different models for the 
second (8542A) CaT line in the left panel and 
the third (8662A) CaT line in the right panel. 
For the models with [Fe/H]< -2, non-LTE ef- 
fects make significant contributions to the line- 
strength. The gray dotted lines indicate a sim- 
ple fit to these values, using the metallicity and 
temperature of the model as input parameters. 
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Fig. 3. The log(g) and effective temperature space occupied 
by a selection of the atmospheric models (dots) and theoreti- 
cal isochrones (lines). In the left panel Yonsei-Yale isochrones 
are shown for different metallicities ([Fe/H]=-1 dashed lines, 
[Fe/H]=-2 solid lines, and [Fe/H]=-3 dotted lines) and ages (12 
Gyr, 8 Gyr, and 2 Gyr (respectively blue, red and black)). In the 
right panel the Yonsei-Yale isochrones for 12 Gyr and [Fe/H]=-1 
and [Fe/H] =-2 are plotted (solid lines) and overplotted with the 
same isochrones from BaSTI (dashed lines). All isochrones plot- 
ted here assume an alpha-element enhancement of [a'/Fe]=+0.4. 



imum and minimum values for the equivalent widths and that 
these errors are systematic. 

3. Because our synthetic grid is not a stellar system, the 
HB magnitude does not have an obvious meaning. Therefore, 
to compare our models with empirical relations which require 
the height above the HB as an input parameter, we have to rely 
on observational or theoretical relations between My of the HB 
( Vhb) of a system and its metallicity. The My for each of the grid 
RGB stars is taken from the isochrones. This value is consistent 
with the value we get if we calculate My from Mboi u s ing th e 
bolometric correction of giant stars from Alonso et al,| (Il99 9|). 
To det ermine Vhb we use the relation given in Catelan & Cortisl 
(l2008h which is calculated using theoretical models for RR 
Lyrae stars. Within its uncertainties, this relation is in excel- 
lent agreement with ob servations of the Vhb of globular clusters 
(e.g.. (Rich etal]l2005h . 
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In order to make a fair comparison between the synthetic 
spectra and the empirical relation, the atmospheric model prop- 
erties were chosen to be as close as possible to known globu- 
lar cluster properties. In the atmospheric models we therefore 
chose to use the alpha-enhanced models with [a/Fe] =-1-0.4, ex- 
cept at the higher metallicities ([Fe/H]>-1 .0) where the MARCS 
grid only provides lower [a/Fe] to match observations in the 
Milky Way. Furthermore, for all [or/Fe] =-1-0.4 models, [Ca/Fe] 
was set to -1-0.25 in the Turbospectrum program, which even 
more closely resembles the true values obser ved in the Galactic 
halo globular clusters (e.g.. lPritzl et al.ll2005l) . For the isochrone 
set we use an age of 12 Gyr, comparable to measured ages for 
the globular clusters (e.g., Krauss & Chaboyer 2 003.) . We find 
that for a range of old ages, between 8 and 15 Gyr, the exact 
choice of the isochrone age does not significantly affect our re- 
sults. 

The successful comparison between our synthetic spectra 
and the empirical relations, using the BaSTI and the Yonsei-Yale 
sets of isochrones, is shown in Fig. |4] The results for the two 
different sets of isochrones are comparable, although the BaSTI 
isochrones give a slightly better coverage and agreement at the 
tip of the RGB. In general, we obtain a good match between the 
predictions of our synthetic spectra equivalent widths and the 
empirical relation, especially for the most luminous part of the 
RGB and at intermediate metallicities where the relation is best 
studied. At [Fe/H]=-0.5 the match is clearly worse, but also the 
empirical relation is not well constrained at this metallicity. We 
also find a larger deviation from t he empirical relati on closer to 
the HB. It was alreadv predicted bv lPont et"al] (|2004|) using mod- 
els that the strength of the CaT lines incre ases more rapidly for 
the more luminous part of the RGB. Later. ICarrera et all (l2007h 
reported a nonlinear tendency in the equivalent width versus 
absolute magnitude relation at fainter magnitudes from an ob- 
servational study of a large sample of RGB stars in open and 
globular cluste r s over a wide range of magnitudes. Recently, 
iDa Costa et al.l (l2009h reported a flattening of the slope of the 
relation below the HB in two globular clusters. These observa- 
tions and early predictions appear to confirm the effect we see in 
our synthetic spectra. This trend is observed to be even stronger 
below the HB (Carrera et al. 200?.; Da Costa et al.. ,2009) . and 
clearly shows one has to be very cautious applying any of the 
[Fe/H]-CaT relations to faint stars, especially below the HB. In 
this paper, we only focus on the RGB above the HB. 



4.2. Further calibration 

In addition to the comparison with the existing empirical re- 
lations, we also calibrate our spectral synthesis models with 
two (very) well studied examples, namely the Sun using the 
Kurucz solar flux atlas (Kurucz et al. 1984) and Arcturus us- 
ing the Hinkle Arcturus atlas ( Hinkle et aL 2000). Although the 
match for the Sun is very good, the initial comparison between 
the observational Hinkle Arcturus spectrum and the synthesized 
spectrum from our models was not satisfactory, especially in the 
wings of the strongest CaT line (A 8542 A). Because this line 
is extremely broad, it is possible that some of the outer parts of 
the line were mistakenly taken to be the continuum level dur- 
ing the continuum subtraction. After careful renormalization of 
the continuum at this wavelength region we were able to get 
an acceptable match using the abundances for Arcturus from 
iFulbright et aH (l2007h . 
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Fig. 5. From synthetic analysis we measure the equivalent widths 
of the two strongest CaT lines as in Fig. |4] extended to lower 
metallicities. Color coding and symbols for -2.5 < [Fe/H] < 
-0.5 is the same as in Fig. |4] the additional metallicities are 
[Fe/H]=-2.5 (purple squares), [Fe/H]=-3.0 (blue plus signs), 
and [Fe/H]=^.0 (pink crosses). The empirical relations (black 
dash-dotted lines), including Icr uncertainties (gray dash-dotted 
lines) are shown for [Fe/H] =-0.5, [Fe/H] =-1.0, [Fe/H] =-1.5, 
[Fe/H]=-2.0, [Fe/H]=-2.5, and [Fe/H]=-3.0 (from top to bot- 
tom). 

5. The CaT lines at [Fe/H] < -2.5 

5.1. The empirical relation 

Given the success in reproducing the well established calibra- 
tion of CaT in the range -2.0 < [Fe/H] < -0.5, we now extend 
our synthetic spectral analysis down to [Fe/H]=^.0. The results 
are shown in Fig.|5]together with the empirical relation extended 
linearly to the low-metallicity regime. We use the Yonsei-Yale 
isochrone set since it extends down to the lowest metallicitiefl 
From a comparison of the empirical relation (dash-dotted lines) 
and the synthetic spectra (colored symbols) in Fig. |5j it is ob- 
vious that the match with the empirical linear relations breaks 
down (as expected) at low metallicities, starting from [Fe/H]<- 
2.0. For example, the empirical relation for [Fe/H]=-3.0 lies be- 
low the synthetic spectra predictions for [Fe/H] =-4.0 models for 
fainter RGB stars. 

As can be seen in Fig.|5] there is also a trend with the inferred 
height of a 'star' above the HB (e.g., surface gravity and temper- 
ature) and the deviation of the 'observed' metallicity compared 

' The Yonsei-Yale isochrone set provides isochrones with [Fe/H]>- 
3.6, so we have to use the [Fe/H]=-3.6 isochrone for our [Fe/H]=^ 
models. Because at lower metallicities the isochrones come closer to- 
gether the error introduced by this mismatch is negligible. 
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Fig. 4. We plot equivalent widths of the two strongest CaT lines which are measured in the synthetic spectra (symbols) and the 
empirical relation from Battaglia et al. (2008b) (black dash-dotted lines) versus (V-Vhb), calculated as described in the text. For 
the empirical relation we derived Icr uncertainties from the spread of the individual RGB stars in globular clusters from Battaglia 
et al. (2008b) around the empirical relation, these are shown as gray dash-dotted lines. The synthetic spectra are calculated for 
metallicities [Fe/H]=-0.5 (gray asterisks), [Fe/H]=-1.0 (black diamonds), [Fe/H]=-1.5 (red triangles) and [Fe/H]=-2.0 (green 
circles). The empirical relation is calculated for the same range of metallicities (from top to bottom). The synthetic spectra are 
matched to the physical properties of RGB stars using the BaSTI isochrones (left panel) and the Yonsei-Yale isochrones (right 
panel). The error bars denote the uncertainties due to the calculation of HB magnitudes (horizontally) and the log(g)-Teff relation 
from isochrones (vertically), for which the error bars show a maximum and minimum value. 



to the metallicity of the model, for the low-metallicity models. 
This demonstrates the fact that not just the equivalent width, but 
also the slope of the relation is changing. At lower metallicities, 
the equivalent width of the line becomes less sensitive to varia- 
tions in gravity or temperature of the star (i.e., its position on the 
RGB). 

The mismatch between the extended empirical relation and 
the synthetic modeling predictions at low metallicities is empha- 
sized in Fig. |6] where the input metallicity of the models is plot- 
ted versus the metallicity obtained from the empirical relation as 
given in Eqs.[T]and|2l From Fig.|6] we obtain valuable insight 
into how extremely low-metallicity spectra would appear in, for 
instance, the DART sample of classical dwarf galaxies. While 
some of the models with input [Fe/H] < -2.5 are correctly re- 
produced by the synthetic spectral method, there are also exam- 
ples where the metallicity is seriously overestimated. The linear 
empirical relation thus offers no means to discriminate between 
very low or extremely low metallicity RGB stars. 

6. A new calibration 

To re-calibrate the relation between the CaT equivalent width 
and metallicity for the low-metallicity regime, we use the equiv- 
alent widths ob tained from ou r synthe tic spectra. In the rela- 
tion obtained bv lBattaglia et alj (l2008bh . there are three free pa- 
rameters to fit the slope of the relation as a function of height 



above the HB, expressed in W , and a linear relation between 
W and [Fe/H]. To extend this calibration we expect to need at 
least two more parameters to fit the dominant features of the 
low-metallicity models in one relation: the changing slope with 
[Fe/H] and the changing offset between lines of equal metal- 
licity. As input into our fitting routine we use the results from 
the synthetic spectra grid from models described in Table [T] with 
[ff/Fe] = -hO.4 (and [Ca/Fe] set to -1-0.25). We use the Yonsei-Yale 
set of isochrones, because only those provide the low metallici- 
ties needed for the a nalysis. For the fitting we use the IDL func- 
tion MPFIT2DFUN dMarkw ardt 2009) to fit a plane relating the 
absolute magnitude and equivalent widths of the synthetic mod- 
els to their metallicity and an F-test to distinguish the best fit. 
Not all models are given the same weights in the fit, for the 
models both higher up the RGB ((V - Vhb)<-1) and at interme- 
diate metallicities (-1.5<[Fe/H]<-0.5) we give higher weights 
since these models represent the best studied region of, and 
best matching models with, the empirical relation. Additionally, 
all models at the (extremely) low-metallicity regime we are 
particularly interested in (below [Fe/H] =-2) are also weighted 
more. On the other hand, the highest metallicity models (at 
[Fe/H] =-0.5) are weighted less, since these have higher uncer- 
tainties and fit less well the empirical relation. Several functional 
forms have been explored. The best fit was obtained by a sim- 
ple linear relation in both luminosity and equivalent width, as in 
the classical empirical calibrations, with merely two extra terms. 
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Fig. 6. From synthetic analysis, we plot the 'real' input value for 
[Fe/H] in the model versus the 'observed' value of [Fe/H] for the 
grid of models, obtained by treating the model spectra as if they 
were observed in the DART program. For each 'real' [Fe/H] in- 
put value are several points representing RGB stars of the same 
metallicity at several places along the RGB. Error bars are calcu- 
lated from the uncertainties shown in Fig.|4]and described in the 
text. If the CaT method would work perfectly all models should 
fall (within their uncertainties) on the one-to-one relation shown 
by the solid black line. Clearly there is an increasing deviation 
starting at [Fe/H]<-2. Same symbols and color coding used as 
in Fig. m 

One term of equivalent width to the power -1.5 to account for 
the variations at low metallicities, one cross term to account for 
the slight non-linear slope, as given in Eq.|5j 



[Fe/H] = -2.87 + 0.195 x (V - Vhb) + 0.458 x EW(2+3) 

-0.913 X EW(2'4) 4- 0.0155 x EW(2+3) x (V - Vhb) (5) 

This relation is only calibrated for RGB stars above the HB 
and should thus not be applied to stars outside -3 < (V - Vhb) < 
0. We want to stress that this relation remains 'empirical', in the 
sense that no theoretical arguments are used to find the best fit- 
ting formula. The corresponding fit is shown in Fig.|2l it fits both 
the higher metallicity end (and in this regime the existing empir- 
ical relation) and the lower metallicity models well within their 
uncertainties. From Fig. [8] which shows the residual of the fit in 
the right panel, it can be seen that the new calibration performs 
less well at the high-metallicity end of the calibration ([Fe/H]^- 
0.5), but even there the error is still within a typical observational 
error bar for [Fe/H]. We estimate the typical maximum error on 
the fitted parameters in Eq.|5]to be ~8%, on the basis of Monte- 
Carlo simulations of the uncertainties on (V-Vhb) and the equiv- 
alent widths. These reasonably low error values convince us that 
the parameters in our new CaT calibration are quite robust to 
changes in our approximations. 

Additionally, in Appendix |A] we describe the relation be- 
tween [Fe/H], EW and My and Mi, to enable the use of the CaT 
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Fig. 7. Same as Fig.|5] but with our new calibration overplotted 
as thick solid black lines. 



[Fe/H]., 




Fig. 8. Two density plots showing the [Fe/H] values for the com- 
plete set of models (left panel) and the difference between the 
[Fe/H] values of the models and their calculated [Fe/H] using 
the new calibration (right panel). 



lines as a metallicity estimator for individual RGB stars in sys- 
tems without a well-defined horizontal branch or for individual 
RGB field stars. 



6.1. Verifying tine new caiibration at iow-metaiiicity 

To verify the reliability of our models at very low-metallicities 
([Fe/H]<-2.0), we have measured the CaT equivalent widths 
for six low-metallicity halo stars with existing high-resolution 
spectroscopic analyses. The properties of these stars are given 
in Table [31 Their CaT spectra are degraded to a resolution of 
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Table 2. The parameters of the observed stars used for verification of the modeling results at low metallicities 



9 



star 


CaT reference 


HR reference 


Teff (K) 


log(g) (cm/s^) 


[Fe/H] HR (dex) 


[Fe/H] CaT(new) (dex) 


HD122563 


W. Aoki, priv. comm. 


Honda et al. (2004) 


4570 


1.1 


-2.77 


-2.66 


HD110184 


J. Fulbright, priv. comm. 


Fulbright (2000) 


4400 


0.6 


-2.3 


-2.27 


HD88609 


Marrese et al. f2003') 


Fulbright (2000) 


4450 


0.6 


-2.9 


-2.97 


HD4306 


J. Fulbright, priv. comm. 


Fulbright (2000) 


4800 


1.7 


-2.8 


-2.57 


HD216143 


J. Fulbright, priv. comm. 


Fulbright (2000) 


4525 


1.0 


-2.1 


-2.10 


CD -38 245 


VLT X-shooter commissioning 


Cavrel et al. (2004) 


4800 


1.5 


-4.2 


-3.82 


Bool 1137 


Norris et al. (2008) 


Norris et al. (2008) 
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Fig, 9. Plotted with our new CaT calibration (colored lines) are 
the well-studied RGB stars in the Milky Way halo (black dia- 
monds) and one extremely low-metallicity star in Bootes I (black 
asterisk) which are all described in Table |2] 



R=6500, equal to the resolution of the synthetic spectra grid, 
and their CaT EWs are measured as described in Section 14.11 
My for these halo stars is calculated from Mboi using the spec- 
troscopically defined values for gravity an d temperature and th e 
bolometric correction of giant stars from lAlonso et al.l (Il999h . 
Additionally their height above the horizontal branch can be ap- 
proximated using their spectroscopic metallicity and the relation 
from Catelan & Cortes (2008). In Fig.|9]the results of applying 
our new CaT calibration to these observations are shown. There 
is clearly a very good agreement. 

Additi onally, one s t ar from the Bootes I dwarf galaxy as 
studied by iNorris et all (l2008l l2009h using medium- and high- 
resolution spectroscopy is also plotted. For this star only one of 
the t wo strongest CaT lines could be measured by iNorris et aH 
(1200 8) with confidence, and the total equivalent width for both 
lines was inferred using t his single line and the obser ved ratio be- 
tween the two fines from lNorris et al.l (12008 . 1996^1 Norris et al.1 
( |2008 |) find [Fe/H] =-3 .45 from medium-resolution spectroscopy 
using the Ca II H and K lines, which is very close to the value 
we deduce from the EW of the CaT line in the same spec- 
trum, [Fe/H]=-3.32. In their subsequent high-resolution follow- 
up study [Fe/ H] is measured dir ectly from Fe lines, which gives 
[Fe/H]=-3.7 (INorris et al.ll2009h . 



6.2. The DART low-metallicity follow-up program 

As a complementary approach to determine if there are any 
extremely low metallicity stars in the dwarf galaxies, DART 
has undertaken a follow-up program to obtain HR spectroscopy 
using the Subaru Telescope High Dispersion Spectograph, 
the UVES spectrograph at VLT, and the MIKE spectrograph 
at Magellan for a sample of stars with CaT [Fe/H] j -2.5 
(lAoki et al.1 l2009l Tafelmeyer et al. in prep., Venn et al. in 
prep.). These HR spectra were taken as an addition to the al- 
ready existing HR spectra from the main program of DART us- 
ing VL T/FLAMES with the GIRAFFE spect ograph in Medusa 
mode dAoki et al.1 120091: iTolstov et all l2009l Letarte et al. in 
press. Hill et al. in prep., Venn et al. in prep.). In these follow-up 
programmes several extremely low-met allicity stars have been 
found, with [Fe/H] values below -3.0 dAoki et al .1120091 Venn 
et al. in prep., Tafelmeyer et al. in prep.) and even around ^.0 
(Tafelmeyer et al. in prep.). Figure [TO] shows all the HR results, 
compared to their LR [Fe/H] values inferred from the CaT lines 
using both the old and the new calibration. In this figure the lim- 
iting range of the old (linear) calibration is also clearly visible, 
only the new calibration extends down to the lowest metallici- 
ties. For the lower metallicities the error bars become larger, due 
to the fact that the relations for different metallicities lie closer 
together and thus a similar error in equivalent width results in 
a larger error in [Fe/H]. Taking this into account, the new cali- 
bration appears to give an accurate prediction of the HR [Fe/H] 
values. 

There are however some stars showing a deviation larger 
than Icr, where the most clear example is the extremely metal- 
poor star in Fornax. This deviation might be (partly) due to non- 
LTE effects or other deficiencies in the modeling of the HR spec- 
trum in order to derive the stellar parameters and abundances for 
this really low-gravity star Some support of this explanation is 
the fact that the agreement between the LR (CaT) and HR re- 
sults improves when [Fell/H] is used instead of [Fel/H], most 
clearly for the extremely metal-poor Fornax star (Tafelmeyer 
et al. in prep.). Non-LTE eff'ects are gener ally negligible for 
the dominant ionization state of Fe II (e.g., Thevenin & Idiarg 



ll999HKraft & Ivansl2003l:lMashonkina et al.l2009h . however for 
Fe I the non-LTE effects are expected to be more significant in 



low-m etallicity, lo w-gravity stars (e.g ., Thevenin & Idiartll 19991; 



^ From our synthetic modeling results we can check this value by 
evaluating the ratio in line strength between the two strongest CaT lines. 



These results are shown in Appendix IB] 



iGehren et al..,2001clMashonkinaet alj|2009h 



7. Alpha element dependence on the CaT lines 

It is naturally expected that differences in [Ca/Fe] will signifi- 
cantly change the equivalent width of the CaT lines, and there- 
fore alter the observed relation between equivalent width and 
[Fe/H]. The MARCS collaboration also provides model atmo- 
spheres where the a elements are not enhanced, but are kept sim- 
ilar to the ratio in the Sun for all models with [Fe/H] > -2. Figure 
[TT]shows that we recover this difference in our synthetic spectra. 
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Fig. 10. A comparison between the CaT [Fe/H] and the [Fe/H] values as determined from high-resolution measurements for the 
same stars from the DART survey, for both the old (left panel) and new CaT calibrations (right panel). Stars from the galaxies 
Sculptor (black circles), Fornax (red asterisks). Carina (blue diamonds) and Sextans (green triangles) are included. 



Qualitatively, the space between the enhanced and solar [Ca/Fe] 
models at equal [Fe/H] seems to agree with the step taken in 
abundance. The equivalent width in our grid of synthetic spec- 
tra does trace the Ca abundance. Since the strength of any line 
depends on the line opacity divided by the continuous opacity, it 
is expected that also other elements than Ca can affect the CaT 
line EW. If these elements contribute free electrons these can en- 
hance the concentration and therefore affect the continuous 
opacity. Which element contributes most free electrons is depen- 
dent on the effective temperature and layer of the atmosphere of 
the star in consideration, but for cool stellar atmospheres in gen- 
eral tiiejnain_sources_of electrons are Mg, Fe, Si, Ca, Na, and 
Al (IShetrone et alJl2009 b. Since thus both Fe and some of the 
a elements have to be considered as important electron donors, 
this can significantly affect the dependence of the CaT EW on 
the [a/Fe] ratio. Nevertheles we find that, within the abundance 
parameters we adopt in this study, the Ca abundance itself is by 
far the dominant factor determining the EW of the CaT line, as 
can be seen in Figure (TT] 

However, when comparing the sensitivity of the CaT equiv- 
alent width measurements of RGB stars in Sculptor and 
Fornax with h i gh-reso lution [Fe/H] and [Ca/H] measurements 
iBattagUa etaP fcoOSbh find that the CaT equivalent width is ac- 
tually a more robust estimator of Fe than Ca. This would suggest 
that it is therefore not advisa ble to u se the CaT as a linear esti- 
mator for [Ca/H] (see Bat taglia et al.ll2008 b. and their Figs. 10 
and 12). This result is not expected from the theoretical expec- 
tations shown in Fig.[TT] There are several factors that may con- 
tribute to this apparent discrepancy between the modeling and 
observational results. First, in our grid of synthetic spectra we 
assume the relative abundance of Ca to the other a elements to 
stay constant, an assumption that not necessarily holds for all 



RGB stars. The extra free electrons donated by the other a ele- 
ments might affect the strength of th e CaT line through th e con- 
tinuous opacity as described above. I Shetrone et all 0009 ) find 
that a 0.5 dex decrease in the electron contributors can increase 
the CaT line strength enough to mimic a Ca abundance increase 
of ~0.2 or even ~0.4 depending on the atmosphere parameters 
of the star. However, in Fornax we generally find that the im- 
portant electron donors Fe and Mg are increased relatively to 
the values in our grid at similar [Ca/H], while we still measure 
the CaT lines to be too strong relative to the [Ca/H] measure- 
ments directly from Cal lines. This result clearly indicates that 
the effect of electron donors can not be driving the offset be- 
tween [Ca/H] derived from CaT and HR Cal analyses. Second, 
the HR Ca abundances are usually derived from fewer Ca 1 lines 
compared to the large number of HR Fe 1 lines available and are 
therefore subject to larger observational errors. Third, the HR 
deter mination of the Cal ab undance is subject to non-LTE ef- 
fects (Mashonkina et al. 2007). If non-LTE effects are included 
in the analysis of the HR spectra to determine the Ca abundances 
this might lead to a closer match in Ca abundances derived from 
the Cal and CaT lines. To investigate this more closely we have 
modeled both the Ca 1 and CaT line strengths using abundances 
and atmospheric parameters from the well-studied halo star CD 
-38 245 (e.g., Cavrel et al. 2004) as a test case. The modeling, 
using the same models and techniques as described in Sect. 13.21 
is performed for [Ca/Fe] =+0.4 in LTE and non-LTE to deter- 
mine the offset from non-LTE to LTE abundances. The results 
are shown in Table [3] It can be seen that the LTE approximation 
has an effect on the determination of the Ca abundance from the 
Ca I lines (a difference of over 0.2 dex for the Ca I line at /I = 
6162 A). In the LTE approximation, the agreement between the 
abundances derived from the Cal and CaT lines is very poor. 
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Fig. 11. Synthetic spectra model predictions for [Ca/Fe]=+0.25 
and [Ca/Fe]=0.0. Models with the same [Fe/H] are plotted us- 
ing the same symbols and colors. However, the synthetic spec- 
tra models with [Ca/Fe]=+0.25 are plotted with larger sym- 
bols and are connected by dotted lines whereas the models with 
[Ca/Fe]=0.0 have smaller symbols and are connected by dashed 
lines. The [Fe/H] and [Ca/H] values for each set of models are 
indicated. 



Table 3. Ca abundances for CD -38 245 



element 


wavelength (A) 


A(Ca) LTE 


A(Ca) non-LTE 


Cal 


4226.73 


2.17 


2.15 


Cal 


6162.18 


2.12 


2.35 


Call 


8498.02 


3.97 


3.10 


Call 


8542.09 


3.42 


2.83 


Call 


8662.14 


3.69 


2.92 
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Fig. 12. Equivalent widths vs. the height above the HB for the 
DART dataset members of Sculptor, Fornax, Carina and Sextans 
(gray circles) in separate panels and the new CaT calibration at 
different [Fe/H] (colored lines). Typical error bars for the ob- 
served DART RGB stars for three luminosity bins are given in 
the upper part of each panel. 



CaT analysis results in a Ca abundance much greater than for 
Ca I. In non-LTE, we are able to reproduce better all the mea- 
sured Ca features in an extremely metal-poor star. The remaining 
discrepancy probably relates to the outer atmospheric layers that 
are not very well modeled - even in non-LTE. To fully resolve 
the discrepancy between the Cal and CaT results, one would 
have to properly explore the effects of specific details including 
line profile fitting and uncertainties in the stellar parameters. 

The stars in the dataset described bv lBattaglia et al.l (l2008bh 
which show the largest discrepancies in CaT compared to 
[Ca/H], are typically much more metal -rich RGB stars than CD - 
38 245. For this regime we have not modeled the non-LTE effects 
in Cal lines. Most of the discrepant stars are from the Fornax 
dwarf galaxy, where we were just able to target its brightest pop- 
ulation of RGB stars due to the relatively large distance to this 
dwarf galaxy (see also Section ing . This means we are also sta- 
tistically probing closer to the tip of the RGB, where we find the 
lower gravity stars for which the outer layers are more diffuse 
and thus more difficult to model - certainly assuming LTE. 



8. Implications for the DART survey 

Our new CaT- [Fe/H] calibration enables a more direct search 
for extremely metal-poor stars in existing datasets, like the 
large DART sample of CaT measurements in RGB stars of 
four classical dwarf galaxies (Sculptor, Fornax, Carina, and 
Sextans). These DART samples were observed in LR (R~ 6500) 
Medusa mode using the European Southern Observatory (ESO) 
V LT/FLAMES facility (Pasauini et al. 200j) an d are described 
in 'Helmi et all (12006^) . ,Koch et al., (.20061) and iBattaglia etaH 
(2006, 2008a.y). 

In Fig.[T2]we overplot the complete low-resolution samples 
for the four galaxies observed in the DART program. Sculptor, 
Fornax, Carina and Sextans, along with the new CaT calibration 
(colored lines). All observed RGB stars, which are likely to be 
members, from DART are shown as small gray circles. All stars 
have a signal-to-noise ratio larger than 10 (per A), a velocity er- 
ror smaller than 5 km/s and a velocity which is within 3cr of the 
systemic radial velocity of the galaxy (for Fornax 2.5cr is used, 
because of the larger contribution from Milky Way foreground^ 
These criteria are identical to those applied by Battaglia et al] 
(2006, 2008a ). We use Vhb for the dwarf galaxies as given 
in Jrwin & Hatzidimitrioul (Il995h . For the errors in the sum of 
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equivalent widths we use the results of iBattaglia et al ] (l2Q08bh 
who found that the random error from repeated measurements of 
the low-resolution sample in the sum of the two broadest lines 
is well represented by crEW2^3~ 6/(S/N). Although this error 
can be quite extensive (mean error bars per galaxy and lumi- 
nosity bin are shown in Fig. fTSl i. there are clearly a number of 
stars in these galaxies that are predicted to have [Fe/H]<-3. Our 
low-metallicity candidates need to be followed up with high- 
resolution spectroscopy to verify our prediction. A small sample 
of low-metallicity candidates has already been followed up with 
high-resolution spectrosco py within DART , these results are dis- 
cussed in Sect. l6.2l and in lAoki et al.l (12009). Tafelmeyeret al. (in 
prep.) and Venn et al. (in prep.). The extremely metal-poor stars 
discovered within these DART follow-up programs in Sculptor, 
Fornax, Carina, and Sextans and the one extremely metal-poor 
Scul ptor star discovered by Kirbv et al. (2009) and Frebel et al. 
f2009b already provide an independent confirmation of our pre- 
dictions of the existence of extremely low-metallicity stars in 
these classical dwarf galaxies. 

8.1. Old and new calibration: A comparison 

It should be noted that the total number of stars affected by 
the new calibration in the DART CaT surveys is very small. 
In the current DART data set, just 2.5% of the target stars 
have a metallicity below [Fe/H]=-2.5, using the calibration of 
[Battaglia et al. (2008b). Using the new calibration the fraction 
of low-metallicity stars increases somewhat to about 7.5%. Still 
it is important to realize that ~92.5% of the stars in these systems 
are at a metallicity which is consistent with both the calibration 
of Battaglia et al. ('2008b) and that given in Sect. |6] of this pa- 
per. This implies that the number of low-metallicity candidates 
remains very small compared to the total number of observed 
stars. This result is illustrated in Fig. [T3] where we compare the 
metallicity distributions for each galaxy using two calibrations 
of the C aT data: The calibration for dwarf galaxies published 
by BattagUa et al.l ( t2008b ) (in black), and the re-calibration pre- 
sented here (in red). It is clear that the overall distributions of 
metallicities within these galaxies do not change significantly 
for the different CaT calibrations, the only clear differences are 
in the low-metallicity tails which have become more populated 
and more extended with the new CaT calibration. 



8.2. The low-metallicity tails 

Although their number is small relative to the total number 
of observed stars, the (extremely) metal-poor stars in each of 
these galaxies represents an important and interesting compo- 
nent of the overall stellar population. The low-metallicity pop- 
ulation can reveal details of the chemical evolution of a partic- 
ular galaxy and, by comparison, of the differences and similar- 
ities of the early evolutionary stages of galaxies. In Fig. fT4l we 
show a comparison of the low-metallicity tails from individual 
stars, without application of binning or error estimates, in the 
four classical dwarf galaxies of DART using our new calibra- 
tion for the CaT equivalent widths and in an ensemble of ultra- 
faint dwarf galaxies using low-resolution (R~6500) spectra over 
a large wavelength range compared with an extensive spectral 
library dkirby et al. 2008) and the HES survey of the Milky Way 
halo (ISchorck et al.l2009 ) measuring [Fe/H] from a combination 
of indices for the Ca II K and Ha lines. Since the number of stars 
present, and observed, in these systems can vary quite a lot, a 
crucial step in comparing the low-metallicity tails is to normalize 




0.05 
0.00 



Fig. 13. Distributions of [Fe/H] for the four dwarf spheroidal 
galaxies in the DART survey, using both the old calibration (Eq. 
\T\ and |2] black) and the new calibration (Eq. |5] red) of the CaT 
lines. The error bars are Poissonian. 
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Fig. 14. The low-metallicity tails of the four classical dwarf 
galaxies studied by DART as obtained with the new calibration 
presented here. In the top panel these are normalized at [Fe/H]=- 
2.5 and compared to the low-metallicity tails from the HES sur- 
vey (both corrected and raw) from Schorcketal. (2009) and 
a compilation of ultra-faint dwarf galaxies (Kirbv et al. 20080. 
The bottom panel displays the same low-metallicity tails for the 
classical- and ultra-faint dwarf galaxies, but normalized to the 
total number of observed stars in each system. 
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them at a certain metallicity or to a certain total number of mea- 
surements. In the upper panel of Fig. [14] all the low-metallicity 
tails have been normalized to 1 at [Fe/H] - -2.5. In this way, one 
can compare the shape of the metallicity tails below this value of 
[Fe/H] assuming all surveys are equally complete at this level. In 
the lower panel we have normalized to the total number of stars 
observed in each system to give a feeling for the relative num- 
ber of (extremely) low-metallicity stars observed. Ideally, one 
would like to be able to isolate the first generation of stars in all 
these galaxies and compare these populations. With our current 
understanding it is not possible to make such a clear distinction 
between samples, which is why we resort to the methods shown 
in Fig. [141 However, with future larger samples of stars with ac- 
curate abundances in these galaxies one could envision making 
use of the position of the knee in the [a/Fe] ratios for example 
dTolstov et alJ2009h to select the first generations of stars in each 
galaxy. 

Before interpreting these metallicity tails further we stress 
that there are several caveats which prevent us from a detailed 
analysis of these results at face value. 

First, the errors on the individual measurements are large, 
especially at low [Fe/H] and/or for fainter stars (typical error 
bars are shown in Fig.[T2li and asymmetric. The asymmetry and 
dependence on metallicity arise from the fact that the relation 
between [Fe/H] and equivalent width is not linear at low metal- 
licities, since the lines of equal metallicity are closer together 
A symmetric error bar in the measured equivalent width, there- 
fore results in a much larger error downwards than upwards in 
[Fe/H]. For instance, the one star in the Sextans dwarf galaxy 
that falls below the [Fe/H]=-4 calibration (see Fig.fTSI) will be as- 
signed a very low metallicity using the calibration, which causes 
the Sextans curve in Fig.[T4]to stay above zero at [Fe/H]=^.5. 
However, the upper 1 cr error in equivalent width allows a metal- 
licity of [Fe/H]=-3.8. 

Second, there are several selection and/or sampling effects 
that are difficult to correct for For example, the fact that not all 
stars have been observed in all of these systems. For instance, a 
much larger percentage of the Sextans stars have been observed 
than Fornax stars, although the total sample in Fornax is larger 
(it is a larger galaxy). Also, the absolute magnitude range of the 
populations targeted may not be comparable. From Fig. [12] it 
can clearly be seen that in Fornax, for instance, only stars well 
within 2 magnitudes from the tip of the RGB are observed due 
to its larger distance, while in the other galaxies the samples 
extend to fainter stars. Any luminosity bias that might remain 
in our new calibration would therefore result in a different bias 
for different galaxies in the extent and population of the low- 
metallicity tails. Even if one would observe the same range of 
magnitudes from the tip of the RGB for all galaxies, this could 
still introduce a bias, because of age and metallicity effects on 
the RGB for an extended star formation history. Additionally, 
some galaxies possess metallicity gradients (e.g., iTolstov et alj 
[2004; Battaglia et al. 2006), which means that the final metalUc- 
ity distribution can also be greatly affected by the spatial sam- 
pling within a galaxy. 

Thus, even though from Fig. [T4| one might be tempted to 
conclude that the low-metallicity tail of Fornax is different in 
both shape and number from the other classical dwarf galaxies 
and the halo and it must therefore have had a different chemical 
history in the earliest epochs, it is not at all clear whether such 
a conclusion is justified. It is very well possible that the low- 
metallicity tail of Fornax lies hidden under its very dominant, 
relatively young, metal-rich population and we have not been 



able to retrieve it because we sample a very small fraction of the 
total number of RGB stars in this galaxy. 

However, from the top panel of Fig. [T4| it is clear that the 
significant difference betw een the metal - poo r tail in the halo and 
in the dwarf galaxies froml Helmi et al.l (12006) has disappeared. 
This also brings the metallicity distributi on function for the clas- 
sical dSphs in line with the results of Kirbv et al.l (1200 8) for 
the ultra-faint dwarf galaxies. The shapes of the low-metallicity 
tails, normaUzed at [Fe/H] =-2.5, are now much more similar. 

9. Conclusions 

It is important for our understanding of galactic chemical evolu- 
tion to know whether there are just a few, and in this case how 
many, extremely low-metallicity ([Fe/H]< -3) stars in the clas- 
sical dSph galaxies, or none at all. This can change our view of 
their early evolution and their subsequent role in galaxy forma- 
tion and evolution. Subsequently, this will influence our view on 
the present-day satellite galaxies as possible templates for the 
building blocks which formed the Galactic halo and their rela- 
tion to the more metal-poor ultra-faint dwarf galaxies. 

Using a grid of synthetic spectra based on MARCS atmo- 
sphere models we have investigated the behavior of the three 
CaT absorption lines as metallicity indicators over a range of 
metallicities from [Fe/H] =-0.5 down to [Fe/H]=^. Our models 
agree with the well-known observed linear relations at [Fe/H]>- 
2.0, although small deviations from linearity are fo und closer 
to the HB as were already predicted and observed (iPont et al.l 
200l ICarrera et all l2007h . For [Fe/H]<-2.0, the relations of 
equal metallicity get closer together and are flatter as a func- 
tion of absolute magnitude. This behavior was expected and is 
reflected in the change in the profile of the CaT lines, from wing- 
dominated to core-dominated as the metallicity drops. We have 
provided a new calibration taking this behavior into account, 
which we have successfully tested on several well-known low- 
metallicity giant stars. We have also investigated the role of the 
Ca abundance itself on these lines and found that although the 
lines ought to trace Ca instead of Fe, this often is not observed, 
probably due to inaccuracies in the determination of the HR Ca 
abundance. 

The new calibration of the CaT- [Fe/H] relation described 
here therefore provides a reliable updated method to search for 
low-metallicity stars using just the equivalent widths of two CaT 
lines. This method requires only a small wavelength range of 
the spectrum at low resolution and the absolute magnitude of the 
star 

Applying the new caUbration to the DART data sets shows 
that, although low in relative number, there are several extremely 
low-metallicity RGB candidates in these classical dwarf galaxies 
that deserve further study. This new calibration thus opens pos- 
sibilities to study the (extremely) metal-poor tail that makes up 
only a small, but interesting, fraction of the dwarf galaxy stellar 
populations. 
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Table A.l. Best fitting parameters for the new calibrations 



New Calibration (M„) 



New Colibration (M,) 



Calibrated for [Fe/H] 


Calibrated for [Ca/H] 


Abs. mag. 


parameter 


value 


Abs. mag. 


parameter 


value 


(V-Vhb) 


a 


-2.87 


(V-Vhb) 


a 


-2.62 




b 


0.195 




b 


0.195 




c 


0.458 




c 


0.457 




d 


-0.913 




d 


-0.908 




e 


0.0155 




e 


0.0146 


My 


a 


-2.90 


My 


a 


-2.65 




b 


0.187 




b 


0.185 




c 


0.422 




c 


0.422 




d 


-0.882 




d 


-0.876 




e 


0.0133 




e 


0.0137 


Ml 


a 


-2.78 


Ml 


a 


-2.53 




b 


0.193 




b 


0.193 




c 


0.442 




c 


0.439 




d 


-0.834 




d 


-0.825 




e 


0.0017 




e 


0.0013 



Appendix A: Absolute magnitude [Fe/IH] and [Ca/H] 
calibrations 

In numerous studies the equivalent widths of the CaT lines are 
studied in relation to the height of the stars above the HB. While 
this measure is very convenient for the study of stars in (galactic) 
globular clusters, one might also want to study systems which 
have no well-defined HB magnitude, or even no HB at all. In 
the main body of the paper we have chosen to define an ex- 
tended calibration of the CaT equivalent width as a function of 
(V - Vhb), which also enabled a direct co mparison with the em- 
pirical relation of iBattaglia et al. I (l2008bh . Here we also provide 
a calibration directly to (Johnson-Cousins) My and M/. Note that 
this also means that we can remove the assumption regarding the 
HB magnitude of our modeled 'RGB stars'. An extra advantage 
of using M/ (instead of My) is that it is much less sensitive to age 
effects, and therefore provides more accu rate metallicities for , 
in particular, younger populations of stars jCarrera et al. 
We also provide calibrations for [Ca/H] instead of [Fe/H], using 
our full modeling grid of Table [T] including the models with so- 
lar [a/Fe] values. All calibrations use the same functional form 
given in Eq. lA.ll The parameters for both [Fe/H] and [Ca/H] 
using either (V - Vhb), My, or Mi are given in Table I A. ll 



[Fe/H] or [Ca/H] = a -H b x (Abs. mag.) -1- c x EW(2+3) 

+d X EW^j'^) + e X EW(2+3) x (Abs. mag.) (A. 1) 

Since also these relations are only calibrated for RGB stars 
and above the HB, they should thus not be applied to stars out- 
side -3 < (V - Vhb) < 0, -3 < My < 0.8, or -4 < Mi < 0. 
In Fig. lA.ll we show the results of the calibrations using My 
and Mi for the high-resolution DART sample for both [Fe/H] 
and [Ca/H]. Distance moduli for Sculptor, Fornax, Car i na an d 
Sextans are taken from'Kaluznv et al.' ('l995h. lRizzi etalJ (l2007h . 
[Mateo et al. ( 1998), and Mateo et al. ( 1995) respectively. The re- 
sults do not change significantly depending on which absolute 
magnitude is used, as can also be seen from a comparison of 
the upper panels of Fig. lA. II with Fig.[TOl It remains clear from 
Fig. lA.ll that our new CaT calibrations based on the synthesized 
spectra trace the HR [Fe/H] values much better than HR [Ca/H], 
as was already discuss ed in Sect. Q of this paper and found by 
iBattagha et all (l2008bh . 




Fig. A.l. Top panels: Same as Fig.[TO] but now on the y-axis the 
new calibration for My and Mi are plotted. Bottom panels; The 
new caUbration for [Ca/H] using either My or Mi. 



Appendix B: Ratio of the 8542A to the 8662A CaT 
line 

Another useful application of our models is the investigation of 
the ratios between the line strengths of the CaT lines, as a func- 
tion of [Fe/H] and luminosity. This enables the determination 
of the summed equivalent width if just one of the two strongest 
lines are present or usable in the spectrum of a star Additionally, 
investigation of the relative differences in development of the 
lines for various combinations of atmospheric parameters can 
give valuable insights into the line formation process. In Fig. 
IB. II the ratios between the EWs of the second and third CaT 
lines, both of which are first corrected for non-LTE effects using 
Eqs. [3]and|4l are plotted as a function of [Fe/H] for a synthetic 
'giant' spectrum. The variation in the ratio with absolute magni- 
tude turns out to be very small as can be seen from the dispersion 
of the equal metallicity models. Although for the intermediate 
metallicity (-3<[Fe/H]<-l) spectra the ratio seems roughly con- 
stant with metallicity, the ratio changes significantly at higher 
([Fe/H]=-0.5), or extremely low ([Fe/H]=^) metallicities. We 
also find that the value for intermediate metallicities is some- 
what lower than the value of 1.34 as measured by Norris et al] 
(l2008h . Our mean ratio from the models is 1 .27 if all models - 
2.5<[Fe/H]<-0.5 are taken into account (this particular range is 
chosen since it represents the metallicities of the vast majority 
of stars in the classical dwarf spheroidal galaxies). This value 
is confirmed in a comparison with the observed spectra from the 
DART dataset where a very similar mean value is found although 
the spread is much larger due to the observational uncertainties. 
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Fig.B.l. The ratio of the second (8542A) to the third (8662A) 
CaT line as a function of [Fe/H]. The color coding and symbols 
for the different [Fe/H] values are the same as in Fig. |5] 



References 

Alonso, A., Arribas, S., & Martmez-Roger, C. 1999, A&AS, 140, 261 
Alvarez, R. & Plez, B. 1998, A&A, 330, 1109 
Aoki, W., Arimoto, N., Sadakane, K., et al. 2009, A&A, 502, 569 
Ai-mandroff, T. E. & Da Costa, G. S. 1991, AJ, 101, 1329 
Armandroff, T. E. & Zinn, R. 1988, AJ, 96, 92 

Barklem, R S., Christlieb, N., Beers, T. C, et al. 2005, A&A, 439, 129 
Battaglia, G., Helmi, A., Tolstoy, E., et al. 2008a, ApJ, 681, L13 
Battaglia, G., Irwin, M., Tolstoy, E., et al. 2008b, MNRAS, 383, 183 
Battaglia, G., Tolstoy, E., Helmi, A., et al. 2006, A&A, 459, 423 
Butler, K. & Giddings, J. R. 1985, Newsletter of Analysis of Astronomical 
Spectra 9 

Carrera, R., Gallart, C., Pancino, E., & Zinn, R. 2007, AJ, 134, 1298 

Castelli, F. & Kumcz, R. L. 2003, in lAU Syinposium, Vol. 210, Modelling of 

Stellar Atmospheres, ed. N. Piskunov, W. W. Weiss, & D. F. Gray, 20P-+ 
Catelan, M. & Cortes, C. 2008, ApJ, 676, L135 
Cayrel, R., Depagne, E., Spite, M., et al. 2004, A&A, 416, 1 1 17 
Cenarro, A. J., Cardiel, N., Gorgas, J., et al. 2001, MNRAS, 326, 959 
Cenarro, A. J., Gorgas, J., Cardiel, N., Vazdekis, A., & Peletier, R. F 2002, 

MNRAS, 329, 863 
Cohen, J. G. 1978, ApJ, 221, 788 
Cohen, J. G. & Huang, W. 2009, ApJ, 701, 1053 

Cole, A. A., Smecker-Hane, T. A., Tolstoy, E., Bosler, T. L., & Gallagher, J. S. 

2004, MNRAS, 347, 367 
Da Costa, G. S., Held, E. V., Saviane, I., & Gullieuszik, M. 2009, ApJ, 705, 1481 
Demarque, P, Woo, J., Kim, Y, & Yi, S. K. 2004, ApJS, 155, 667 
Frebel, A., Kirby, E., & Simon, J. D. 2009, ArXiv e-prints 
Frebel, A., Simon, J. D., Geha, M., & Willman, B. 2010, ApJ, 708, 560 
Freeman, K. & Bland-Hawthorn, J. 2002, ARA&A, 40, 487 
Fulbright, J. P 2000, AJ, 120, 1841 

Fulbright, J. P, McWilliam, A., & Rich, R. M. 2007, ApJ, 661, 1152 

Fulbright, J. P, Rich, R. M., & Castro, S. 2004, ApJ, 612, 447 

Gehren, T, Bufler, K., Mashonkina, L., Reetz, J., & Shi, J. 2001, A&A, 366, 981 

Giddings, J. R. 1981, PhD thesis, AA(Moscow state University) 

Gustafsson, B., Edvardsson, B., Eriksson, K., et al. 2008, A&A, 486, 951 

Helmi, A., Irwin, M. J., Tolstoy, E., et al. 2006, ApJ, 651, L121 

Hinkle, K., Wallace, L., Valenti, J., & Harmer, D. 2000, Visible and Near Infrared 

Adas of the Arcturus Spectrum 3727-9300 A (ASP) 
Honda, S., Aoki, W., Kajino, T, et al. 2004, ApJ, 607, 474 
Irwin, M. & Hatzidimitriou, D. 1995, MNRAS, 277, 1354 
Ivans, 1. 1., Kraft, R. P, Sneden, C, et al. 2001, AJ, 122, 1438 
Jones, J. E., Alloin, D. M., & Jones, B. J. T. 1984, ApJ, 283, 457 



Jorgensen, U. G., Carlsson, M., & Johnson, H. R. 1992, A&A, 254, 258 
Kaluzny, J., Kubiak, M., Szymanski, M., et al. 1995, A&AS, 112, 407 
Kirby, E. N., Guhathakurta, P, Bolte, M., Sneden, C, & Geha, M. C. 2009, ApJ, 
705, 328 

Kirby, E. N., Simon, J. D., Geha, M., Guhathakurta, P, & Frebel, A. 2008, ApJ, 
685, L43 

Koch, A., Grebel, E. K., Gilmore, G. R, et al. 2008a, AJ, 135, 1580 
Koch, A., Grebel, E. K., Kleyna, J. T, et al. 2007b, AJ, 133, 270 
Koch, A., Grebel, E. K., Wyse, R. F G., et al. 2006, AJ, 131, 895 
Koch, A., McWilliam, A., Grebel, E. K., Zucker, D. B., & Belokurov, V. 2008b, 
ApJ, 688, L13 

Koch, A., Wilkinson, M. 1., Kleyna, J. T, et al. 2007a, ApJ, 657, 241 
Koch, A., Wilkinson, M. 1., Kleyna, J. T., et al. 2009, ApJ, 690, 453 
Kraft, R. P & Ivans, 1. 1. 2003, PASP 115, 143 
Krauss, L. M. & Chaboyer, B. 2003, Science, 299, 65 

Kupka, F. G., Ryabchikova, T. A., Piskunov, N. E., Stempels, H. C, & Weiss, 

W. W. 2000, Baltic Astronomy, 9, 590 
Kurucz, R. 1993, ATLAS9 Stellar Atmosphere Programs and 2 km/s 

grid. Kurucz CD-ROM No. 13. Cambridge, Mass.: Smithsonian 

Astrophysical Observatory, 1993., 13 
Kurucz, R. L., Furenlid, I., Brault, J., & Testerman, L. 1984, Solar flux atlas from 

296 to 1300 nm (National Solar Observatory Atlas, Sunspot, New Mexico: 

National Solar Observatory) 
Leaman, R., Cole, A. A., Venn, K. A., et al. 2009, ApJ, 699, 1 
Markwardt, C. B. 2009, in Astronomical Society of the Pacific Conference 

Series, Vol. 411, Astronomical Society of the Pacific Conference Series, ed. 

D. A. Bohlender, D. Durand, & P Dowler, 251-+ 
Marrese, P M., Boschi, R, & Munari, U. 2003, A&A, 406, 995 
Mashonkina, L., Gehren, T, Shi, J., Kom, A., & Grupp, P. 2009, ArXiv e-prints 
Mashonkina, L., Kom, A. J., & Przybilla, N. 2007, A&A, 461, 261 
Mateo, M., Rischer, P, & Krzeminski, W. 1995, AJ, 1 10, 2166 
Mateo, M., Hurley-Keller, D., & Nemec, J. 1998, AJ, 115, 1856 
Munoz, R. R., Majewski, S. R., Zaggia, S., et al. 2006, ApJ, 649, 201 
Norris, J. E., Rreeman, K. C, & Mighell, K. J. 1996, ApJ, 462, 241 
Norris, J. E., Gilmore, G., Wyse, R. R G., et al. 2008, ApJ, 689, LI 13 
Norris, J. E., Yong, D., Gilmore, G., & Wyse, R. F G. 2009, Ai'Xiv e-prints 
Olszewski, E. W., Schommer, R. A., Suntzeft', N. B., & Harris, H. C. 1991, AJ, 

101, 515 

Pasquini, L., Avila, G., Blecha, A., et al. 2002, The Messenger, 110, 1 
Pietrinferni, A., Cassisi, S., Salaris, M., & Castelli, F 2004, ApJ, 612, 168 
Plez, B. 1998, A&A, 337, 495 

Plez, B. 2008, Physica Scripta Volume T, 133, 014003 
Pont, F, Zinn, R., Gallart, C, Hardy, E., & Winnick, R. 2004, AJ, 127, 840 
Pritzl, B. J., Venn, K. A., & Irwin, M. 2005, AJ, 130, 2140 
Rich, R. M., Corsi, C. E., Cacciari, C, et al. 2005, AJ, 129, 2670 
Rizzi, L., Held, E. V, Saviane, 1., Tully, R. B., & Gullieuszik, M. 2007, MNRAS, 
380, 1255 

Rutledge, G. A., Hesser, J. E., & Stetson, P B. 1997a, PASP 109, 907 
Rutledge, G. A., Hesser, J. E., Stetson, P B., et al. 1997b, PASP 109, 883 
Sbordone, L., Bonifacio, P., Castelli, R, & Kurucz, R. L. 2004, Memorie della 

Societa Astronomica Italiana Supplement, 5, 93 
Schorck, T, Christlieb, N., Cohen, J. G., et al. 2009, A&A, 507, 817 
Shetrone, M., Venn, K. A., Tolstoy, E., et al. 2003, AJ, 125, 684 
Shetrone, M. D., Cote, P, & Sargent, W. L. W. 2001, ApJ, 548, 592 
Shetrone, M. D., Siegel, M. H., Cook, D. O., & Bosler, T. 2009, AJ, 137, 62 
Simon, J. D. & Geha, M. 2007, ApJ, 670, 313 
Spinrad, H. & Taylor, B. J. 1969, ApJ, 157, 1279 
Spinrad, H. & Taylor, B. J. 1971, ApJS, 22, 445 
Spite, M., Cayrel, R., Plez, B., et al. 2005, A&A, 430, 655 
Suntzefl', N. B., Mateo, M., Terndrup, D. M., et al. 1993, ApJ, 418, 208 
Thevenin, F & Idiart, T. P 1999, ApJ, 521, 753 
Tolstoy, E., Hill, V, Irwin, M., et al. 2006, The Messenger, 123, 33 
Tolstoy, E., Hill, V, & Tosi, M. 2009, ARA&A, 47, 371 
Tolstoy, E., Irwin, M. J., Helmi, A., et al. 2004, ApJ, 617, LI 19 
Tolstoy, E., Venn, K. A., Shetrone, M., et al. 2003, AJ, 125, 707 
Venn, K. A., Irwin, M., Shetrone, M. D., et al. 2004, AJ, 128, 1177 
Walker, M. G., Belokurov, V, Evans, N. W., et al. 2009b, ApJ, 694, L144 
Walker, M. G., Mateo, M., Olszewski, E. W., Sen, B., & Woodroofe, M. 2009a, 

AJ, 137, 3109 

Yi, S., Demarque, P, Kim, Y, et al. 2001, ApJS, 136, 417 



